Logic programming-based Minimal Cut Sets reveal consortium-level therapeutic targets for chronic wound infections.

Minimal Cut Sets (MCSs) identify sets of reactions which, when removed from a metabolic network, disable certain cellular functions. The traditional search for MCSs within genome-scale metabolic models (GSMMs) targets cellular growth, identifies reaction sets resulting in a lethal phenotype if disrupted, and retrieves a list of corresponding gene, mRNA, or enzyme targets. Using the dual link between MCSs and Elementary Flux Modes (EFMs), our logic programming-based tool aspefm was able to compute MCSs of any size from GSMMs in acceptable run times. The tool demonstrated better performance when computing large-sized MCSs than the mixed-integer linear programming methods. We applied the new MCSs methodology to a medically-relevant consortium model of two cross-feeding bacteria, Staphylococcus aureus and Pseudomonas aeruginosa. aspefm constraints were used to bias the computation of MCSs toward exchanged metabolites that could complement lethal phenotypes in individual species. We found that interspecies metabolite exchanges could play an essential role in rescuing single-species growth, for instance inosine could complement lethal reaction knock-outs in the purine synthesis, glycolysis, and pentose phosphate pathways of both bacteria. Finally, MCSs were used to derive a list of promising enzyme targets for consortium-level therapeutic applications that cannot be circumvented via interspecies metabolite exchange.

Development and applications of genome-scale metabolic network models.

The genome-scale metabolic network model is an effective tool for characterizing the gene-protein-response relationship in the entire metabolic pathway of an organism. By combining various algorithms, the genome-scale metabolic network model can effectively simulate the influence of a specific environment on the physiological state of cells, optimize the culture conditions of strains, and predict the targets of genetic modification to achieve targeted modification of strains. In this review, we summarize the whole process of model building, sort out the various tools that may be involved in the model building process, and explain the role of various algorithms in model analysis. In addition, we also summarized the application of GSMM in network characteristics, cell phenotypes, metabolic engineering, etc. Finally, we discuss the current challenges facing GSMM.

The Fur-like regulatory protein MAP3773c modulates key metabolic pathways in Mycobacterium avium subsp. paratuberculosis under in-vitro iron starvation.

Johne's disease (JD) is a chronic enteric infection of dairy cattle worldwide. Mycobacterium avium subsp. paratuberculosis (MAP), the causative agent of JD, is fastidious often requiring eight to sixteen weeks to produce colonies in culture-a major hurdle in the diagnosis and therefore in implementation of optimal JD control measures. A significant gap in knowledge is the comprehensive understanding of the metabolic networks deployed by MAP to regulate iron both in-vitro and in-vivo. The genome of MAP carries MAP3773c, a putative metal regulator, which is absent in all other mycobacteria. The role of MAP3773c in intracellular iron regulation is poorly understood. In the current study, a field isolate (K-10) and an in-frame MAP3773c deletion mutant (ΔMAP3773c) derived from K-10, were exposed to iron starvation for 5, 30, 60, and 90 min and RNA-Seq was performed. A comparison of transcriptional profiles between K-10 and ΔMAP3773c showed 425 differentially expressed genes (DEGs) at 30 min time post-iron restriction. Functional analysis of DEGs in ΔMAP3773c revealed that pantothenate (Pan) biosynthesis, polysaccharide biosynthesis and sugar metabolism genes were downregulated at 30 min post-iron starvation whereas ATP-binding cassette (ABC) type metal transporters, putative siderophore biosynthesis, PPE and PE family genes were upregulated. Pathway analysis revealed that the MAP3773c knockout has an impairment in Pan and Coenzyme A (CoA) biosynthesis pathways suggesting that the absence of those pathways likely affect overall metabolic processes and cellular functions, which have consequences on MAP survival and pathogenesis.

A Machine Learning Approach for Predicting Essentiality of Metabolic Genes.

The identification of essential genes is a key challenge in systems and synthetic biology, particularly for engineering metabolic pathways that convert feedstocks into valuable products. Assessment of gene essentiality at a genome scale requires large and costly growth assays of knockout strains. Here we describe a strategy to predict the essentiality of metabolic genes using binary classification algorithms. The approach combines elements from genome-scale metabolic models, directed graphs, and machine learning into a predictive model that can be trained on small knockout data. We demonstrate the efficacy of this approach using the most complete metabolic model of Escherichia coli and various machine learning algorithms for binary classification.

Stable States of a Microbial Community Are Formed by Dynamic Metabolic Networks with Members Functioning to Achieve Both Robustness and Plasticity.

A more detailed understanding of the mechanisms underlying the formation of microbial communities is essential for the efficient management of microbial ecosystems. The stable states of microbial communities are commonly perceived as static and, thus, have not been extensively examined. The present study investigated stabilizing mechanisms, minority functions, and the reliability of quantitative ana-lyses, emphasizing a metabolic network perspective. A bacterial community, formed by batch transferred cultures supplied with phenol as the sole carbon and energy source and paddy soil as the inoculum, was analyzed using a principal coordinate ana-lysis (PCoA), mathematical models, and quantitative parameters defined as growth activity, community-changing activity, community-forming activity, vulnerable force, and resilience force depending on changes in the abundance of operational taxonomic units (OTUs) using 16S rRNA gene amplicon sequences. PCoA showed succession states until the 3rd transferred cultures and stable states from the 5th to 10th transferred cultures. Quantitative parameters indicated that the bacterial community was dynamic irrespective of the succession and stable states. Three activities fluctuated under stable states. Vulnerable and resilience forces were detected under the succession and stable states, respectively. Mathematical models indicated the construction of metabolic networks, suggesting the stabilizing mechanism of the community structure. Thirteen OTUs coexisted during stable states, and were recognized as core OTUs consisting of majorities, middle-class, and minorities. The abundance of the middle-class changed, whereas that of the others did not, which indicated that core OTUs maintained metabolic networks. Some extremely low abundance OTUs were consistently exchanged, suggesting a role for scavengers. These results indicate that stable states were formed by dynamic metabolic networks with members functioning to achieve robustness and plasticity.

Using design of experiments to guide genetic optimization of engineered metabolic pathways.

Design of experiments (DoE) is a term used to describe the application of statistical approaches to interrogate the impact of many variables on the performance of a multivariate system. It is commonly used for process optimization in fields such as chemical engineering and material science. Recent advances in the ability to quantitatively control the expression of genes in biological systems open up the possibility to apply DoE for genetic optimization. In this review targeted to genetic and metabolic engineers, we introduce several approaches in DoE at a high level and describe instances wherein these were applied to interrogate or optimize engineered genetic systems. We discuss the challenges of applying DoE and propose strategies to mitigate these challenges. ONE-SENTENCE SUMMARY: This is a review of literature related to applying Design of Experiments for genetic optimization.

Transposon sequencing reveals metabolic pathways essential for Mycobacterium tuberculosis infection.

New drugs are needed to shorten and simplify treatment of tuberculosis caused by Mycobacterium tuberculosis. Metabolic pathways that M. tuberculosis requires for growth or survival during infection represent potential targets for anti-tubercular drug development. Genes and metabolic pathways essential for M. tuberculosis growth in standard laboratory culture conditions have been defined by genome-wide genetic screens. However, whether M. tuberculosis requires these essential genes during infection has not been comprehensively explored because mutant strains cannot be generated using standard methods. Here we show that M. tuberculosis requires the phenylalanine (Phe) and de novo purine and thiamine biosynthetic pathways for mammalian infection. We used a defined collection of M. tuberculosis transposon (Tn) mutants in essential genes, which we generated using a custom nutrient-rich medium, and transposon sequencing (Tn-seq) to identify multiple central metabolic pathways required for fitness in a mouse infection model. We confirmed by individual retesting and complementation that mutations in pheA (Phe biosynthesis) or purF (purine and thiamine biosynthesis) cause death of M. tuberculosis in the absence of nutrient supplementation in vitro and strong attenuation in infected mice. Our findings show that Tn-seq with defined Tn mutant pools can be used to identify M. tuberculosis genes required during mouse lung infection. Our results also demonstrate that M. tuberculosis requires Phe and purine/thiamine biosynthesis for survival in the host, implicating these metabolic pathways as prime targets for the development of new antibiotics to combat tuberculosis.

Genome-wide characterization of circulating metabolic biomarkers.

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8-11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.

Machine learning predicts system-wide metabolic flux control in cyanobacteria.

Metabolic fluxes and their control mechanisms are fundamental in cellular metabolism, offering insights for the study of biological systems and biotechnological applications. However, quantitative and predictive understanding of controlling biochemical reactions in microbial cell factories, especially at the system level, is limited. In this work, we present ARCTICA, a computational framework that integrates constraint-based modelling with machine learning tools to address this challenge. Using the model cyanobacterium Synechocystis sp. PCC 6803 as chassis, we demonstrate that ARCTICA effectively simulates global-scale metabolic flux control. Key findings are that (i) the photosynthetic bioproduction is mainly governed by enzymes within the Calvin-Benson-Bassham (CBB) cycle, rather than by those involve in the biosynthesis of the end-product, (ii) the catalytic capacity of the CBB cycle limits the photosynthetic activity and downstream pathways and (iii) ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) is a major, but not the most, limiting step within the CBB cycle. Predicted metabolic reactions qualitatively align with prior experimental observations, validating our modelling approach. ARCTICA serves as a valuable pipeline for understanding cellular physiology and predicting rate-limiting steps in genome-scale metabolic networks, and thus provides guidance for bioengineering of cyanobacteria.

Prediction and integration of metabolite-protein interactions with genome-scale metabolic models.

Metabolites, as small molecules, can act not only as substrates to enzymes, but also as effectors of activity of proteins with different functions, thereby affecting various cellular processes. While several experimental techniques have started to catalogue the metabolite-protein interactions (MPIs) present in different cellular contexts, characterizing the functional relevance of MPIs remains a challenging problem. Computational approaches from the constrained-based modeling framework allow for predicting MPIs and integrating their effects in the in silico analysis of metabolic and physiological phenotypes, like cell growth. Here, we provide a classification of all existing constraint-based approaches that predict and integrate MPIs using genome-scale metabolic networks as input. In addition, we benchmark the performance of the approaches to predict MPIs in a comparative study using different features extracted from the model structure and predicted metabolic phenotypes with the state-of-the-art metabolic networks of Escherichia coli and Saccharomyces cerevisiae. Lastly, we provide an outlook for future, feasible directions to expand the consideration of MPIs in constraint-based modeling approaches with wide biotechnological applications.

Bacteria use a catabolic patchwork pathway of apparently recent origin for degradation of the synthetic buffer compound TRIS.

The synthetic buffer compound TRIS (2-amino-2-(hydroxymethyl)propane-1,3-diol) is used in countless applications, and no detailed information on its degradation has been published so far. Herein, we describe the discovery of a complete bacterial degradation pathway for TRIS. By serendipity, a Pseudomonas strain was isolated from sewage sludge that was able to grow with TRIS as only carbon and nitrogen source. Genome and transcriptome analyses revealed two adjacent gene clusters embedded in a mobile genetic element on a conjugative plasmid to be involved in TRIS degradation. Heterologous gene expression revealed cluster I to encode a TRIS uptake protein, a TRIS alcohol dehydrogenase, and a TRIS aldehyde dehydrogenase, catalyzing the oxidation of TRIS into 2-hydroxymethylserine. Gene cluster II encodes a methylserine hydroxymethyltransferase (mSHMT) and a d-serine dehydratase that plausibly catalyze the conversion of 2-hydroxymethylserine into pyruvate. Conjugational plasmid transfer into Pseudomonas putida KT2440 enabled this strain to grow with TRIS and with 2-hydromethylserine, demonstrating that the complete TRIS degradation pathway can be transmitted by horizontal gene transfer. Subsequent enrichments from wastewater purification systems led to the isolation of further TRIS-degrading bacteria from the Pseudomonas and Shinella genera carrying highly similar TRIS degradation gene clusters. Our data indicate that TRIS degradation evolved recently via gene recruitment and enzyme adaptation from multiple independent metabolic pathways, and database searches suggest that the TRIS degradation pathway is now globally distributed. Overall, our study illustrates how engineered environments can enhance the emergence of new microbial metabolic pathways in short evolutionary time scales.

Unraveling the role of the mitochondrial one-carbon pathway in undifferentiated thyroid cancer by multi-omics analyses.

The role of the serine/glycine metabolic pathway (SGP) has recently been demonstrated in tumors; however, the pathological relevance of the SGP in thyroid cancer remains unexplored. Here, we perform metabolomic profiling of 17 tumor-normal pairs; bulk transcriptomics of 263 normal thyroid, 348 papillary, and 21 undifferentiated thyroid cancer samples; and single-cell transcriptomes from 15 cases, showing the impact of mitochondrial one-carbon metabolism in thyroid tumors. High expression of serine hydroxymethyltransferase-2 (SHMT2) and methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) is associated with low thyroid differentiation scores and poor clinical features. A subpopulation of tumor cells with high mitochondrial one-carbon pathway activity is observed in the single-cell dataset. SHMT2 inhibition significantly compromises mitochondrial respiration and decreases cell proliferation and tumor size in vitro and in vivo. Collectively, our results highlight the importance of the mitochondrial one-carbon pathway in undifferentiated thyroid cancer and suggest that SHMT2 is a potent therapeutic target.

BioISO: An Objective-Oriented Application for Assisting the Curation of Genome-Scale Metabolic Models.

As the reconstruction of Genome-Scale Metabolic Models (GEMs) becomes standard practice in systems biology, the number of organisms having at least one metabolic model is peaking at an unprecedented scale. The automation of laborious tasks, such as gap-finding and gap-filling, allowed the development of GEMs for poorly described organisms. However, the quality of these models can be compromised by the automation of several steps, which may lead to erroneous phenotype simulations. Biological networks constraint-based In Silico Optimisation (BioISO) is a computational tool aimed at accelerating the reconstruction of GEMs. This tool facilitates manual curation steps by reducing the large search spaces often met when debugging in silico biological models. BioISO uses a recursive relation-like algorithm and Flux Balance Analysis (FBA) to evaluate and guide debugging of in silico phenotype simulations. The potential of BioISO to guide the debugging of model reconstructions was showcased and compared with the results of two other state-of-the-art gap-filling tools (Meneco and fastGapFill). In this assessment, BioISO is better suited to reducing the search space for errors and gaps in metabolic networks by identifying smaller ratios of dead-end metabolites. Furthermore, BioISO was used as Meneco's gap-finding algorithm to reduce the number of proposed solutions for filling the gaps.

Genome-scale metabolic network models for industrial microorganisms metabolic engineering: Current advances and future prospects.

The construction of high-performance microbial cell factories (MCFs) is the centerpiece of biomanufacturing. However, the complex metabolic regulatory network of microorganisms poses great challenges for the efficient design and construction of MCFs. The genome-scale metabolic network models (GSMs) can systematically simulate the metabolic regulation process of microorganisms in silico, providing effective guidance for the rapid design and construction of MCFs. In this review, we summarized the development status of 16 important industrial microbial GSMs, and further outline the technologies or methods that continuously promote high-quality GSMs construction from five aspects: I) Databases and modeling tools facilitate GSMs reconstruction; II) evolving gap-filling technologies; III) constraint-based model reconstruction; IV) advances in algorithms; and V) developed visualization tools. In addition, we also summarized the applications of GSMs in guiding metabolic engineering from four aspects: I) exploring and explaining metabolic features; II) predicting the effects of genetic perturbations on metabolism; III) predicting the optimal phenotype; IV) guiding cell factories construction in practical experiment. Finally, we discussed the development of GSMs, aiming to provide a reference for efficiently reconstructing GSMs and guiding metabolic engineering.

Flux sampling in genome-scale metabolic modeling of microbial communities.

BACKGROUND: Microbial communities play a crucial role in ecosystem function through metabolic interactions. Genome-scale modeling is a promising method to understand these interactions and identify strategies to optimize the community. Flux balance analysis (FBA) is most often used to predict the flux through all reactions in a genome-scale model; however, the fluxes predicted by FBA depend on a user-defined cellular objective. Flux sampling is an alternative to FBA, as it provides the range of fluxes possible within a microbial community. Furthermore, flux sampling can capture additional heterogeneity across a population, especially when cells exhibit sub-maximal growth rates. RESULTS: In this study, we simulate the metabolism of microbial communities and compare the metabolic characteristics found with FBA and flux sampling. With sampling, we find significant differences in the predicted metabolism, including an increase in cooperative interactions and pathway-specific changes in predicted flux. CONCLUSIONS: Our results suggest the importance of sampling-based approaches to evaluate metabolic interactions. Furthermore, we emphasize the utility of flux sampling in quantitatively studying interactions between cells and organisms.

Genome-scale metabolic modelling of extremophiles and its applications in astrobiological environments.

Metabolic modelling approaches have become the powerful tools in modern biology. These mathematical models are widely used to predict metabolic phenotypes of the organisms or communities of interest, and to identify metabolic targets in metabolic engineering. Apart from a broad range of industrial applications, the possibility of using metabolic modelling in the contexts of astrobiology are poorly explored. In this mini-review, we consolidated the concepts and related applications of applying metabolic modelling in studying organisms in space-related environments, specifically the extremophilic microbes. We recapitulated the current state of the art in metabolic modelling approaches and their advantages in the astrobiological context. Our review encompassed the applications of metabolic modelling in the theoretical investigation of the origin of life within prebiotic environments, as well as the compilation of existing uses of genome-scale metabolic models of extremophiles. Furthermore, we emphasize the current challenges associated with applying this technique in extreme environments, and conclude this review by discussing the potential implementation of metabolic models to explore theoretically optimal metabolic networks under various space conditions. Through this mini-review, our aim is to highlight the potential of metabolic modelling in advancing the study of astrobiology.

Genome-scale metabolic network model and phenome of solvent-tolerant Pseudomonas putida S12.

BACKGROUND: Pseudomonas putida S12 is a gram-negative bacterium renowned for its high tolerance to organic solvents and metabolic versatility, making it attractive for various applications, including bioremediation and the production of aromatic compounds, bioplastics, biofuels, and value-added compounds. However, a metabolic model of S12 has yet to be developed. RESULTS: In this study, we present a comprehensive and highly curated genome-scale metabolic network model of S12 (iSH1474), containing 1,474 genes, 1,436 unique metabolites, and 2,938 metabolic reactions. The model was constructed by leveraging existing metabolic models and conducting comparative analyses of genomes and phenomes. Approximately 2,000 different phenotypes were measured for S12 and its closely related KT2440 strain under various nutritional and environmental conditions. These phenotypic data, combined with the reported experimental data, were used to refine and validate the reconstruction. Model predictions quantitatively agreed well with in vivo flux measurements and the batch cultivation of S12, which demonstrated that iSH1474 accurately represents the metabolic capabilities of S12. Furthermore, the model was simulated to investigate the maximum theoretical metabolic capacity of S12 growing on toxic organic solvents. CONCLUSIONS: iSH1474 represents a significant advancement in our understanding of the cellular metabolism of P. putida S12. The combined results of metabolic simulation and comparative genome and phenome analyses identified the genetic and metabolic determinants of the characteristic phenotypes of S12. This study could accelerate the development of this versatile organism as an efficient cell factory for various biotechnological applications.

GINv2.0: a comprehensive topological network integrating molecular interactions from multiple knowledge bases.

Knowledge bases have been instrumental in advancing biological research, facilitating pathway analysis and data visualization, which are now widely employed in the scientific community. Despite the establishment of several prominent knowledge bases focusing on signaling, metabolic networks, or both, integrating these networks into a unified topological network has proven to be challenging. The intricacy of molecular interactions and the diverse formats employed to store and display them contribute to the complexity of this task. In a prior study, we addressed this challenge by introducing a "meta-pathway" structure that integrated the advantages of the Simple Interaction Format (SIF) while accommodating reaction information. Nevertheless, the earlier Global Integrative Network (GIN) was limited to reliance on KEGG alone. Here, we present GIN version 2.0, which incorporates human molecular interaction data from ten distinct knowledge bases, including KEGG, Reactome, and HumanCyc, among others. We standardized the data structure, gene IDs, and chemical IDs, and conducted a comprehensive analysis of the consistency among the ten knowledge bases before combining all unified interactions into GINv2.0. Utilizing GINv2.0, we investigated the glycolysis process and its regulatory proteins, revealing coordinated regulations on glycolysis and autophagy, particularly under glucose starvation. The expanded scope and enhanced capabilities of GINv2.0 provide a valuable resource for comprehensive systems-level analyses in the field of biological research. GINv2.0 can be accessed at: https://github.com/BIGchix/GINv2.0 .

The Reactome Pathway Knowledgebase 2024.

The Reactome Knowledgebase (https://reactome.org), an Elixir and GCBR core biological data resource, provides manually curated molecular details of a broad range of normal and disease-related biological processes. Processes are annotated as an ordered network of molecular transformations in a single consistent data model. Reactome thus functions both as a digital archive of manually curated human biological processes and as a tool for discovering functional relationships in data such as gene expression profiles or somatic mutation catalogs from tumor cells. Here we review progress towards annotation of the entire human proteome, targeted annotation of disease-causing genetic variants of proteins and of small-molecule drugs in a pathway context, and towards supporting explicit annotation of cell- and tissue-specific pathways. Finally, we briefly discuss issues involved in making Reactome more fully interoperable with other related resources such as the Gene Ontology and maintaining the resulting community resource network.

Chemical-tag-based semi-annotated metabolomics facilitates gene identification and specialized metabolic pathway elucidation in wheat.

The importance of metabolite modification and species-specific metabolic pathways has long been recognized. However, linking the chemical structure of metabolites to gene function in order to explore the genetic and biochemical basis of metabolism has not yet been reported in wheat (Triticum aestivum). Here, we profiled metabolic fragment enrichment in wheat leaves and consequently applied chemical-tag-based semi-annotated metabolomics in a genome-wide association study in accessions of wheat. The studies revealed that all 1,483 quantified metabolites have at least one known functional group whose modification is tailored in an enzyme-catalyzed manner and eventually allows efficient candidate gene mining. A Triticeae crop-specific flavonoid pathway and its underlying metabolic gene cluster were elucidated in further functional studies. Additionally, upon overexpressing the major effect gene of the cluster TraesCS2B01G460000 (TaOMT24), the pathway was reconstructed in rice (Oryza sativa), which lacks this pathway. The reported workflow represents an efficient and unbiased approach for gene mining using forward genetics in hexaploid wheat. The resultant candidate gene list contains vast molecular resources for decoding the genetic architecture of complex traits and identifying valuable breeding targets and will ultimately aid in achieving wheat crop improvement.